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ABSTRACT 

Results are presented from a survey for molecular anions in seven nearby Galactic star-forming cores and molecular 
clouds. The hydrocarbon anion CgH - is detected in all seven target sources, including four sources where no 
anions have been previously detected: LI 172, L1389, L1495B, and TMC-1C. The Cg H - /C f , H column density ratio 
is >1.0% in every source, with a mean value of 3.0% (and standard deviation 0.92%). Combined with previous 
detections, our results show that anions are ubiquitous in dense clouds wherever Cg H is present. The CgH ~/CgH 
ratio is found to show a positive correlation with molecular hydrogen number density, and with the apparent age of 
the cloud. We also report the first detection of C 4 H - in TMC-1 (at 4.8er confidence), and derive an anion-to-neutral 
ratio C 4 H-/C 4 H = (1.2 ± 0.4) x 10 -5 (= 0.0012% ± 0.0004%). Such a low value compared with CgH - highlights 
the need for a revised radiative electron attachment rate for C 4 H. Chemical model calculations show that the 
observed C 4 H could be produced as a result of reactions of oxygen atoms with CgH and CgH - . 
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1. INTRODUCTION 

Negative ions (anions) were first discovered in the interstellar 
medium (ISM) by McCarthy et al. (2006), who detected the 
linear hydrocarbon anion CgH - in TMC-1. This was followed 
by CgH - and C 4 H - detections in the protostellar core L1527 by 
Sakai et al. (2007) and Agundez et al. (2008), respectively, which 
hinted of a role for anions in the chemistry of star formation. 
Gupta et al. (2009) performed the first dedicated survey for 
CgH - in 24 molecular sources, but detected the anion in only two 
star-forming clouds, with anion-to-neutral ratios on the order of 
a few percent, similar to previously observed values. Interstellar 
anion detections were confined to the Taurus molecular cloud 
complex until the recent discoveries of CgH - in the Lupus, 
Cepheus, and Auriga star-forming regions (Sakai et al. 2010; 
Cordiner et al. 2011), which proved that anions are widespread 
in the local ISM, and not an artifact of any particular physical 
or chemical conditions in the Taurus region. 

The possible importance of anions in interstellar chemistry 
was first discussed by Dalgarno & McCray (1973), and prospects 
for the detection of molecular anions using radio astronomy 
were examined by Sarre (1980). Herbst (1981) argued that 
large interstellar molecules (including carbon-chain-bearing 
species) can undergo rapid radiative electron attachment — as 
exemplified by the laboratory experiments of Woodin et al. 
(1980) — potentially resulting in significant anion abundances 
in dense molecular clouds. Based on this idea, models for anion 
chemistry have been successful in reproducing the observed 
abundances of CgH - and C 8 H - in TMC-1, IRC+ 102 16, LI 527, 
and L1512 (Millar et al. 2007; Remijan et al. 2007; Harada & 
Herbst 2008; Cordiner et al. 2008, 2012). However, there is a 
major discrepancy between the modeled and observed C 4 H - 
anion-to-neutral ratio (see Herbst & Osamura 2008), and the 
lack of anions in PDRs (Agundez et al. 2008) is at variance 
with the model predictions of Millar et al. (2007). Clearly, 
our understanding of molecular anion chemistry is incomplete. 


Nevertheless, as demonstrated by Cordiner et al. (2012), anion 
measurements have the potential to offer insight into interstellar 
cloud properties due to their reactivity (Eichelberger et al. 
2007), and consequently, their sensitivity to the abundances of 
gas-phase electrons and atomic C, H, and O, and to physical 
conditions such as cloud density and molecular depletion. 
Observations of anions in a variety of interstellar environments 
will be key to a complete understanding of their chemistry. Given 
the relatively small number of molecular anion detections in the 
ISM to date, we set out to address the question of just how 
widespread anions are in interstellar clouds, and to ascertain the 
behavior of the anion-to-neutral ratio over a range of cloud types 
and ages. 

In this article, results are presented from our survey for CgH - 
and CgH in a sample of seven carbon-chain-rich interstellar 
clouds, protostars, and prestellar cores. Detections of CgH - 
in two of the sources (LI 251 A and L 1 5 12) were previously 
reported by Cordiner et al. (2011), and here we report the 
final results for all seven sources, as well as the first (tentative) 
detection of C 4 H - in TMC-1. 

2. TARGET SOURCES 

Due to the close chemical relationship between polyynes and 
cyanopolyynes (see, e.g., Federman et al. 1990; Millar & Herbst 
1994), carbon-chain-bearing species such as C 4 H and CgH are 
expected to be abundant in dense molecular clouds, close to 
the peak of HC 3 N emission. Target sources with strong HC 3 N 
emission lines were selected from an HC 3 N J = 10-9 mapping 
survey of 20 nearby molecular clouds, prestellar cores, and 
young, low-mass protostars obtained using the Onsala Space 
Observatory (OSO) 20 m telescope between 2005 and 2012. 
Some of these data were presented by Buckle et al. (2006) and 
Cordiner et al. (201 1). Maps of the integrated HC 3 N J — 10-9 
emission intensity toward six of our seven chosen sources 
(LI 172, L1251A, L1389, L1495B, L1512, and TMC-1C) are 
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Table 1 

Targets, Coordinates, and Distances 


Source 

R.A. 

(J2000) 

Decl. 

(J2000) 

Cloud Type a 

Distance 

(pc) 

Ref. 

Ll 172 SMM 

21 : 02 : 22.1 

+67:54:48 

Star forming 

440 

1 

L1251A 

22:30:40.4 

+75:13:46 

Star forming 

300 

2 

L1389 (CB17) SMM1 

04:04:36.6 

+56:56:00 

Protostellar 

250 

3 

L1495B 

04:15:41.8 

+28:47:46 

Quiescent 

140 

4 

L1512 

05:04:07.1 

+32:43:09 

Star forming 

140 

4 

TMC-1C 

04:41:35.6 

+26:00:21 

Star forming 

140 

4 

TMC-1 CP 

04:41:41.9 

+25:41:27 

Quiescent 

140 

4 


Notes. 

a See Section 2 for description of source classification scheme. 

References. (1) Visser et al. 2002; (2) Kun & Prusti 1993; (3) Launhardt et al. 
2010; (4) Elias 1978. 


shown in Figures 1 and 2. These observations were obtained 
using the OSO with a beam size of 42" and beam efficiency of 
0.5 ± 0.05. A 30" map spacing was used, which was resampled 
to a 7" 5 pixel -1 grid using bilinear interpolation. For TMC-1, 
FIC 3 N maps were published by Hirahara et al. (1992) and Pratap 
et al. (1997). 

To search for Cf, H and C’gFT , we targeted the strongest 
HC 3 N peak positions in our maps of L1389, L1495B, L1512, 
and TMC-1C. For LI 172 we targeted the submillimeter core 
cataloged by Di Francesco et al. (2008), which fell within 15" of 
the FIC 3 N peak. Our chosen LI 251 A position does not coincide 
with the main HC 3 N peak because it was based on an earlier, 
lower-sensitivity map than that shown in Figure 1 . For TMC-1, 
we targeted the cyanopolyyne peak (CP) where Agundez et al. 
(2008) previously attempted to detect C 4 H - . The adopted 
coordinates for our anion survey targets are listed in Table 1 . 

Our LI 172 and L1512 Cf,H positions overlap the prestellar 
cores identified by Visser et al. (2002) and Ward-Thompson 
et al. (1994), respectively. For L1389, the HC 3 N emission 
maps out a very compact clump, coincident with the cen- 
ter of the Bok Globule CB17 (Clemens & Barvainis 1988). 
The HC 3 N peak also matches closely the location of strongest 
microwave/ submillimeter continuum emission observed by 
Launhardt et al. (2010) and denoted L1389 SMM1 (the posi- 
tion of which is indicated in the bottom panel of Figure 1 by the 
cross inside the dashed black circle of the Green Bank Telescope 
(GBT) beam). As shown in Figure 1, our targeted LI 25 1 A posi- 
tion is ~40" southeast of the center of the protostar LI 251 IRS 3 
(Lee et al. 2010), although the outer envelope of this protostar 
probably intersects the GBT beam to some extent, as discussed 
by Cordiner et al. (2011). 

From dust continuum emission, Launhardt et al. (2010) ob- 
tained a number density for L1389 SMM1 of 6 x 10 6 cm -3 . 
Pavlyuchenkov et al. (2006) derived a moderately high CO de- 
pletion factor (~40), and inferred a relatively evolved chem- 
ical age of ~2 Myr for this collapsing core. The core was 
recently found to contain a very young, low-luminosity protostar 
(referred to as CB17 MMS; Chen et al. 2012), 8 " NE of our 
targeted GBT beam center, as well a more-evolved Class 0/1 
protostar 20" NW. Given that the radii of typical low-mass pro- 
tostellar envelopes are ~10 4 AU (e.g., Jprgensen et al. 2002), 
at a distance of 250 pc, our targeted LI 389 beam is likely to be 
dominated by protostellar envelope matter. 

L1495B, on the other hand, is a quiescent molecular cloud 
with no known protostars nearby. It appears to be chemically 
young, with an apparently low level of CO depletion consistent 
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Figure 1. Observed HC 3 N J = 10-9 integrated intensity maps toward LI 172, 
L1251A, and L1389. The OSO beam size is indicated in the lower right of 
each panel with a white circle. Intensity scales are given in units of K km s -1 
and have not been corrected for beam efficiency. Dashed black circles represent 
the size (HPBW) and position of our targeted GBT anion survey positions. 
Crosses denote locations of submillimeter (prestellar) cores (Ward-Thompson 
et al. 1994; Visser et al. 2002; Di Francesco et al. 2008; Launhardt et al. 2010) 
and asterisks denote protostellar core locations (Visser et al. 2002; Chen et al. 
2012; Lee et al. 2010). 

(A color version of this figure is available in the online journal.) 
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Figure 2. Same as Figure 1, for L1495B, L1512, and TMC-1C. 
(A color version of this figure is available in the online journal.) 


with less-evolved interstellar clouds (Hirota et al. 2004). In 
TMC-1C, Buckle et al. (2006) identified an anti-correlation 
between the spatial distributions of HC 3 N and C 18 0, at least 
partly attributable to CO depletion. Thus, our observations of 


TMC-1C (at the HC 3 N peak) probably sample more chemically 
evolved, depleted gas than L1495B. No protostars are known 
to be present in our observed TMC-1C beam, but several 
submillimeter sources are located in the surrounding cloud 
(indicated by white crosses in the bottom panel of Figure 2), 
showing this to be a region of active star formation. Southwest 
of TMC-1C lies the source TMC-1 (CP), which is a well-known 
chemically young dark cloud, with only a modest degree of 
depletion (Hirahara et al. 1992; Cordiner et al. 2012). 

In Table 1, we provide a basic categorization of our target 
sources in light of their properties described above: “quiescent” 
refers to those cloud cores that are chemically young, show no 
evidence for active star formation within the telescope beam 
(such as outflows or compact submillimeter/IR emission), and 
do not appear to be undergoing collapse; “star forming” is used 
for clouds showing nearby active star formation and a greater 
degree of chemical evolution; finally, “protostellar” describes 
ourL1389 position, which contains the low-luminosity protostar 
CB17MMS. 


3. ANION OBSERVATIONS 

Observations of emission lines of CgH (/ = 1 0.5-9. 5 ), CgH - 
( J = 10-9 and 11-10), and HC 3 N (J = 3-2 and 4-3) were 
carried out between 2010 April and 2011 January using the Ka 
receiver of the NRAO 100 m GBT 5 . The GBT Spectrometer was 
used with a bandwidth of 50 MHz and 8 1 92 x 6 . 1 kHz channels 
(corresponding to a velocity spacing of ^0.065 km s -1 ). Four 
spectral windows were used, which allowed simultaneous ob- 
servation of the CgH and CgH lines. For the compact, spatially 
isolated source L1512, beam switching (with 78" throw) was 
used, and for all other targets, frequency switching was used. 
Pointing was checked every one to two hours and was typi- 
cally accurate to within 5". In the middle of the observed fre- 
quency range (28 GHz), the telescope beam FWHM was 26" and 
the main beam efficiency was 0.90. Total system temperatures 
were typically in the range 60-80 K and the zenith opacity was 
0.05 ± 0.02. Intensity calibration was performed using beam- 
switched observations of the compact radio source NGC 7027. 
Measured antenna temperatures were subsequently corrected 
for opacity, spillover, ohmic loss, blockage efficiency, and beam 
efficiency, then averaged using standard GBTIDL routines. 

We obtained observations of C 4 H ( N = 2-1, J = 1.5-0. 5, 
F = 2-1) and C 4 H (/ = 2-1) between 2011 December and 
2012 February at around 19 GHz using the GBT A'-band focal- 
plane array. These observations were obtained with channel 
spacing 12.2 kHz, beam efficiency 0.92, beam size 39", and 
zenith opacity 0.025 ±0.015. 

Additional single-pointing HC 3 N J = 10-9 spectra of our 
sources were obtained using the OSO at our adopted target 
positions. A list of all the observed transitions is given in Table 2. 

4. RESULTS 

4.1. New Anion Detections 

The CgH anion and its neutral counterpart CgH were 
detected in all of the surveyed sources. This is the first time 
anions have been detected in LI 172, L1389, L1495B, and 
TMC-1C. The observed CgH - and CgH spectra are shown in 
Figure 3. 


5 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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Table 2 

Observed Species, Transitions, and Frequencies 


Species 

Transition 

Frequency 

(MHz) 

Ref. 

Telescope 

HPBW a 

(") 

c 4 h 

IV = 2-1, J = 1. 5-0.5, F = 2-1 

19054.4760 

1 

GBT 

38 

c 4 h 

J = 2-1 

18619.7580 

2,3 

GBT 

39 

c 6 h- 

J = 10-9 

27537.1302 

4 

GBT 

26 

C 6 H- 

J = 11-10 

30290.8133 

4 

GBT 

24 

c 6 h 

J = 10.5-9.5, /, F = 11-10 

29109.6437 

5 

GBT 

25 

c 6 h 

J = 10.5-9.5, /, F = 10-9 

29109.6855 

5 

GBT 

25 

c 6 h 

J = 10.5-9.5, e, F = 11-10 

29112.7087 

5 

GBT 

25 

c 6 h 

J = 10.5-9.5, e,F = 10-9 

29112.7503 

5 

GBT 

25 

hc 3 n 

J = 3-2, F = 3-3 

27292.9018 

6 

GBT 

27 

hc 3 n 

J = 3-2, F = 2-1 

27294.0758 

6 

GBT 

27 

hc 3 n 

J = 3-2, F = 3-2 

27294.2927 

6 

GBT 

27 

hc 3 n 

J = 3-2, F = 4-3 

27294.3451 

6 

GBT 

27 

hc 3 n 

J = 3-2, F = 2-2 

27296.2334 

6 

GBT 

27 

hc 3 n 

J = 4-3, F = 4-4 

36390.8861 

6 

GBT 

20 

hc 3 n 

J = 4-3, F = 3-2 

36392.2358 

6 

GBT 

20 

hc 3 n 

J = 4-3, F = 4-3 

36392.3293 

6 

GBT 

20 

hc 3 n 

J = 4-3, F = 5-4 

36392.3630 

6 

GBT 

20 

hc 3 n 

J = 4-3, F = 3-3 

36394.1765 

6 

GBT 

20 

hc 3 n 

J = 10-9, F = 9-8 

90978.9838 

6 

OSO 

42 

hc 3 n 

J = 10-9, F = 10-9 

90978.9948 

6 

OSO 

42 

hc 3 n 

J = 10-9, F = 11-10 

90979.0024 

6 

OSO 

42 


Notes. 

a Telescope half-power beam width at observed frequency. 

References. (1) Gottlieb et al. 1983; (2) Gupta et al. 2007; (3) McCarthy & Thaddeus 
2008; (4) McCarthy et al. 2006; (5) McCarthy & Thaddeus 2005; (6) HC 3 N line 
predictions are based on spectroscopic data summarized by Thorwirth et al. (2000) 
and deZafra (1971). 

Spectra of the J = 10-9 and J — 11—10 transitions of C(,H _ 
were averaged in velocity space to improve the signal-to-noise 
ratio; integrated intensities of the averaged spectra are given in 
Table 3. Due to partial blending of the hyperfine components 
of the Q,H (J = 10.5-9.5, /) lines, the integrated intensity 
summed over both components is given in Table 3. Given the 
weakness of these lines, they can safely be assumed to be 
optically thin, and Equation (2) of Lis et al. (2002) was used to 
derive column densities for CgH and C<-,H (given in Table 3). 
Spectroscopic data for the transitions of interest (Table 2) were 
obtained from the Cologne Database for Molecular Astronomy 
(Muller et al. 2005), and an excitation temperature of 10 K was 
assumed. Error estimates were derived using 500 Monte Carlo 
noise replications for each measurement. 

In TMC-1, Q, H excitation temperatures of 5.2 K and 6.7 K 
were measured by Bell et al. (1999) and Sakai et al. (2007), 


respectively. Such sub-thermal excitation might be expected 
given the large Einstein A coefficients for rotational transitions 
of this molecule, which are a direct consequence of its large 
(5.5 D) dipole moment. However, the value obtained by Bell 
et al. (1999) may be subject to uncertainty due to the effects of 
telescope beam dilution, and the Sakai et al. (2007) value may 
not be applicable to our observations because it was obtained 
using the J = 16.5-15.5 transition, which has an Einstein A 
coefficient about four times greater than that of the / = 10.5-9.5 
transition we observed. The larger Einstein A results in an 
increased likelihood of sub-thermal excitation of the J = 16.5 
level. If the C7,H excitation temperature is as low as 5 K, then the 
calculated column densities for C () H and C fl H will be about 
36% smaller than the values in Table 3. Due to the similar 
moments of inertia of Cf,H and Cf,H± their rotational levels 
occur at similar energies, so that the calculated anion-to-neutral 
column density ratios (C^H /C7,H) are insensitive to small 
uncertainties in excitation, provided both molecules share a 
common excitation temperature. The mean C () H ~/Q,H ratio for 
the seven sources is 3.0%, with a standard deviation of 0.92%. 

Our new measurements for C 6 H and Q,H in L1251A and 
LI 5 12 supersede those presented by Cordiner et al. (2011). 
Although the column densities and anion-to-neutral ratios mea- 
sured in the present study match those of Cordiner et al. (201 1) 
within the stated errors, the newer values benefit from improved 
calibration and analysis methods and have smaller uncertain- 
ties. The difference between our TMC-1 (CP) Cf,H~/C () H value 
of 2.5% ±0.4% and the value of 1.6% ±0.3% obtained by 
Briinken et al. (2007) is just beyond the range of the la error 
bars, and could be the result of statistical noise in the spec- 
tra. Alternatively, because Briinken et al. (2007) derived their 
value from a different set of transitions, non-LTE excitation ef- 
fects in CgH and/or C fl H~ could be partially responsible for the 
discrepancy. 

The observed C 4 H - and C 4 H spectra of TMC-1 (CP) are 
shown in the bottom-right panel of Figure 3. The 7 = 2-1 
line of C 4 H- is detected in a single channel with a peak 
antenna temperature of 3.4 mK, which corresponds to 4.8a 
(where a is the rms noise of the baseline-subtracted spec- 
trum). Using a Gaussian fit to this C 4 H line with 1000 Monte 
Carlo noise replications, the central velocity was found to be 
5.70 ± 0.05 km s " 1 with an FWHM of 0.43 ± 0.13 km s" 1 . 
Within the errors, these parameters match those derived for 
HC 3 N (given in Table 4), which adds confidence to our CTH - 
detection. The integrated line intensity is 1 .0 ± 0.3 mK km s -1 , 
which corresponds to a column density of ATC 4 H ") = 
(8.0 ± 2.4) x 10 9 cm' 2 (assuming an excitation temperature 


Table 3 

Molecular Line Measurements and Anion-to-neutral Ratios 


Source 

/ 7M B 4r(C 6 H-) a 
(mK km s -1 ) 

/ 7MB^(C 6 H) b 
(mK km s -1 ) 

N( C 6 H-) 
(10 10 cm" 2 ) 

W(C 6 H) 
(10 11 cm -2 ) 

c 6 h-/c 6 h 

(%) 

Ll 172 

6.7 (0.8) 

41.1 (1.6) 

2.4 (0.3) 

7.1 (0.3) 

3.3 (0.5) 

L1251A 

6.5 (1.0) 

43.6 (2.1) 

2.3 (0.4) 

7.6 (0.4) 

3.0 (0.6) 

L1389 

5.9 (0.8) 

27.1 (1.3) 

2.1 (0.3) 

4.7 (0.2) 

4.4 (0.8) 

L1495B 

9.6 (1.0) 

141.6 (2.2) 

3.4 (0.4) 

24.6 (0.4) 

1.4 (0.2) 

L1512 

4.3 (0.4) 

26.3 (0.9) 

1.5 (0.1) 

4.6 (0.2) 

3.3 (0.4) 

TMC-1C 

13.6(1.1) 

88.1 (1.7) 

4.8 (0.4) 

15.3 (0.3) 

3.1 (0.3) 

TMC-1 CP 

41.6(5.0) 

332.4 (9.0) 

14.7(1.8) 

57.8(1.6) 

2.5 (0.4) 


Notes. Uncertainties given in parentheses are ±68% Monte Carlo errors. Measurements for L1251A and L1512 
supersede those of Cordiner et al. (201 1) as explained in Section 4. 
a Integrated intensities of average C6H - J = 10-9 and J = 11-10 spectra. 

b Integrated intensities summed over both hyperfine components of the C6H J = 10.5-9.5, / transition. 
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Figure 3. Observed CfiH~ spectra (average of J = 10-9 and J = 1 1-10 transitions), and Cr, H J = 10.5-9.5, / spectra. The Q, H velocity scale is given with respect 
to the (weighted) mean frequency of the two hyperfine components. For TMC-1, observed C 4 F 1 - (J = 2-1) and C 4 H (N = 2-1, J = 1.5-0. 5, F = 2-1) spectra are 
also shown (bottom-right panel). Dashed vertical lines indicate the least-squares HC 3 N velocities for each target. 
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Table 4 

HC 3 N RADEX Fit Results 


Source 

/V(HC 3 N) 
(10 13 cm” 2 ) 

u(HC 3 N) 
(km s” 1 ) 

Au(HC 3 N) 
(km s~‘) 

/1 h 2 

(10 4 cm" 3 ) 

Ll 172 

2.7 

2.74 

0.49 

7.5 (2.3) 

L1251A 

3.3 

-3.97 

0.31 

2.1 (0.4) 

L1389 

1.5 

-4.75 

0.32 

5.2 (1.3) 

L1495B 

8.2 

7.62 

0.19 

1.1 (0.3) 

L1512 

4.2 

7.11 

0.10 

2.6 (0.9) 

TMC-1C 

6.4 

5.37 

0.18 

1.1 (0.4) 

TMC-1 CP 

19.5 

5.75 

0.32 

1.0 (0.2) 


Notes. Uncertainties on »h 2 (in parentheses) were calculated assuming ± 1 K 
uncertainty on the kinetic temperature, but may be larger due to the possible 
effects of beam dilution. 


of 10 K). The neutral C 4 H line has an integrated intensity of 
411.3 ± 0.7 mK km s” 1 , which requires A/XC 4 H) = (6.91 ± 
0.01) x 10 14 cm” 2 . The corresponding anion-to-neutral ratio 
C 4 H-/C 4 H is (1.2 ± 0.4) x 10” 5 . This value is consistent with 
the upper limit of5.2x 10 ” 5 measured by Agundezetal. (2008), 
and our C 4 H column density also closely matches their value 
of 7.1 x 10 14 cm” 2 . Both /V(C 4 H) and A/XC 4 H”) are relatively 
insensitive to the adopted excitation temperature, and vary by 
only 12% over the range 5-10 K. 

4.2. HC 3 N Observations and Radiative Transfer Modeling 

Our HC 3 N maps (Figures 1 and 2) show the presence of 
one or more compact molecular condensations in each source, 
and highlight a clumpy structure in this carbon-chain-rich gas. 
The peaks in FIC 3 N emission tend to coincide approximately 
with the locations of known prestellar cores and protostars, 
the positions of which are denoted with crosses and asterisks, 
respectively. As a result of the collapse and infall processes 
occurring in these objects, their densities are expected to be 
greater than the surrounding gas. This should result in increased 
collisional excitation of the 7=10 level of HC 3 N, which has a 
critical density of 7 x 10 5 cm ” 3 (Buckle et al. 2006). Thus, the 
intensity variations shown in our HC 3 N OSO maps represent a 
combination of variations in gas density and total HC 3 N column. 

Spectra of three rotational transitions of HC 3 N for each 
source, observed at the locations targeted by our GBT anion 
survey, are shown in Figure 4. These span a range of upper-state 
energy levels 1.82-16.69 cm ” 1 (2.6-24.0 K). In order to derive 
gas densities and velocities, the HC 3 N spectra were subject to 
fitting using the RADEX radiative transfer code (van der Tak 
et al. 2007). Collisional transition rates were taken from the 
Leiden Atomic and Molecular Database 6 (Schoier et al. 2005), 
which tabulates scaled versions of the original data of Green 
& Chapman (1978). Hyperfine structure (HFS) was accounted 
for by assuming LTE population of the hyperfine levels in each 
J state. The primary collision partner number density («h 2 ), 
molecular column density (/V), Doppler line FWHM (An), and 
Doppler velocity ( v ) were optimized for the spectra of each 
source using the MPFIT least-squares algorithm (Markwardt 
2008). The best-fitting parameters are given in Table 4. The 
kinetic temperature was fixed at 10 K, then varied by ±1 K 
to produce the quoted error estimates on n h 2 . For all sources, 
a temperature of 10 K produced a good fit to all five HFS 
components (the two weakest HFS components are not shown in 
Figure 4). The HC 3 N Doppler velocities are shown with vertical 


6 http://www.strw.leidenuniv.nl/~moldata 


dashed lines in Figure 3 and closely match the central velocities 
of the detected hydrocarbons and anions. In several cases, the 
HC 3 N line widths are very narrow; for L15 12, Au = 0. 1 km s” 1 , 
which is only slightly greater than the 10 K thermal line width 
of 0.07 km s” 1 , indicating an unusual lack of turbulence, flows, 
or shears in the central ~30" of this source. 

The C 6 H”/C 6 H ratio is plotted as a function of n Hl in 
Figure 5, and shows a positive correlation. There is considerable 
scatter and uncertainty on «h 2 , but consistent with previous 
studies (Hirahara et al. 1992; Visser et al. 2002; Launhardt 
et al. 2010; Kirk et al. 2005, respectively), we find the lowest 
density in TMC-1, highest densities in LI 172 and L1389, and 
an intermediate density in L1512. Our calculated n h, values for 
LI 172, L1389, and L1512 are, however, systematically lower 
than previously derived values by up to two orders of magnitude. 
This discrepancy is likely indicative of a problem in our method, 
which can be attributed to the effects of different degrees of 
beam dilution affecting the different HC 3 N line frequencies 
we observed. Given the larger 42" beam size for the HC 3 N 
J = 10-9 line observed with the OSO 20 m compared to the 
27" GBT beam for the J = 3-2 line, it is theoretically possible 
that the OSO data suffer up to an additional factor of 2.6 in beam 
dilution compared with that of the GBT, which could result in 
sufficient reduction of the derived n h 2 values to account for 
the observed discrepancies. Figures 1 and 2 (and Figure 7 of 
Pratap et al. 1997) show that the HC 3 N emission from all our 
observed sources is distributed over an area larger than the OSO 
beam, which suggests that beam dilution may not be so severe. 
However, given the lack of spatial information on the extent of 
the J = 10-9 emission on size scales less than the OSO beam 
size, it is impossible to quantify the importance of beam dilution 
in our RADEX calculations, so the «h 2 values given in Table 4 
should be treated with caution. More accurate estimates for n H2 
could be obtained from HC 3 N J = 10-9 observations using a 
telescope with a smaller beam (such as the IRAM 30 m), which 
would be a better match to that of our GBT observations. 

5. DISCUSSION 
5.1. C 6 H~ 

Since the 2006 discovery of Cr,H” in TMC-1, there have been 
seven reported detections of this anion in various parts of the 
Taurus-Auriga molecular cloud complex (including our latest 
detections in L1495B, L1512, and TMC-1C). Thus far, there 
have been only two detections of interstellar anions outside of 
this small region (Sakai et al. 2010; Cordiner et al. 2011), which 
highlights the importance of our new discoveries of CgH” in 
the vicinity of LI 172 SMM (a prestellar core in Cepheus) and 
LI 389 SMM1 (a prestellar/protostellar core inside Bok Globule 
CB17 in Camelopardalis). These new detections confirm that 
anions are indeed widespread in star-forming regions outside 
of Taurus, and that CeH” appears to be ubiquitous in dense 
interstellar clouds wherever its parent neutral is present. 

Myers et al. (1988) detected an outflow in LI 172 that Visser 
et al. (2002) identified as originating from a protostar (LI 172 
SMM1) that lies about 30" south of our targeted position. They 
also detected a dense core (presumed prestellar), about 50" SE of 
our target, with AThL) = 6 x 10 21 cm” 2 . If this column density 
also applies to our observed position, the C f , H abundance is 
4 x 10” 12 , which is a factor of a few less than in TMC-1, but an 
order of magnitude greater than in L1512, L1521F, and L1544 
(see Cordiner et al. 2012, and references therein), suggesting a 
relatively high absolute anion abundance in LI 172. 
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27.293 27.294 36.391 36.392 36.393 90.978 90.98 


LSR Frequency (GHz) 

Figure 4. Observed HC 3 N spectra (black points) with RADEX least-squares models overlaid (red curves). Frequency-switching residuals have been masked from the 
observed data. Three main hyperfine peaks are visible in the J = 3-2 and J = 4-3 spectra. Model parameters are given in Table 4. 

(A color version of this figure is available in the online journal.) 
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Figure 5. C6H anion-to-neutral ratio vs. H 2 number density ( nu 2 ). The dashed 
line shows linear least-squares fit. 


By utilizing the available physical and chemical information 
on the interstellar clouds for which CgH - has so far been de- 
tected (as summarized for our observed targets in Section 2), 
it is possible to examine the CgH - /CgH ratio with respect to 
each cloud’s chemical and dynamical evolutionary state. The 
evolutionary states of the observed clouds are given in Table 1 , 
and are shown by the different symbols on the plot in Figure 5. 
There is a clear division in the anion-to-neutral ratio between 
young, “quiescent” and older, “star-forming” cores: for the 
quiescent, chemically young cores L1495B and TMC-1 (CP), 
CgH - /CgH < 3%, and for the star-forming cores/protostars, 
CgH - /CgH ^ 3%. L1495B is the most chemically young and 
quiescent of our sources (Hirota et al. 2004), with no evidence 
for nearby star formation, and has the lowest CgH - /CgH ratio in 
our sample (1.4% ± 0.2%). The very low luminosity, low mass 
protostar (VeLLO) LI 389 SMM1, on the other hand, evidently 
contains among the most evolved matter in our sample (having 
passed through the quiescent and star-forming stages), and has 
the highest anion-to-neutral ratio (CgH - /CgH = 4.4% ± 0.8%). 
This value is similar to the value of 4.0% ± 1.0% found in 
L1521F by Gupta et al. (2009), which also contains a VeLLO 
(Bourke et al. 2006). Comparing our results with other sources 
from the literature, L1527 (Sakai et al. 2007) contains a some- 
what more-evolved Class 0/1 protostar, and this source has a 
significantly larger CgH /CgH ratio of 9.3% ± 2.9%. Thus, we 
confirm the apparent relationship identified by Cordiner et al. 
( 201 1 ) between the evolutionary state of interstellar matter and 
its CgH - anion-to-neutral ratio. The quiescent cloud Lupus- 1 A 
observed by Sakai et al. (2010) also fits the trend (with 
CgH - /CgH = 2.1% ± 0.6%), as does the very low lumi- 
nosity protostar Chamaeleon MMS1 (with its upper limit of 
CgH - /CgH < 10%; Cordiner & Charnley 2012). The dense 
prestellar core LI 544 may deviate, however, with a relatively 
low ratio of2.5%±0.8% (Gupta et al. 2009). 

The observed CgH - /CgH trend and its correlation with «h 2 
shown in Figure 5 can be understood in the context of the 
theoretical study of Cordiner et al. (2012), who identified 
the effects of depletion on the CgH - anion-to-neutral ratio, 
the degree of which is related to both the density and chemical 
age of the cloud. Depletion occurs as atoms and molecules 
collide with and stick to dust grains, and thus proceeds faster 
in denser media. Increased depletion affects the anion-to- 
neutral ratio in two ways: ( 1 ) the free electron abundance 


goes up so that more CgH - is produced by radiative electron 
attachment and (2) the atomic O and H abundances go down, 
thus reducing the anion destruction rate. Our CgH - observations 
agree with this theory, particularly for the quiescent cores and 
the protostars, which are believed to lie at opposite ends of 
the interstellar chemical/dynamic evolutionary path. Measured 
anion-to-neutral ratios, however, are not sufficiently accurate 
to draw a clear distinction among the (moderately evolved) 
prestellar cores and star-forming gas clouds of LI 172, L125 1 A, 
L1512, and TMC-1C. These clouds likely contain gas that 
spans a range of densities and degrees of chemical evolution. 
More accurate measurements of density, depletion, and the 
CgH - /CgH ratio at high spatial resolution will be required in 
order to confirm the utility of the anion-to-neutral ratio as a 
measure of the evolutionary state of interstellar clouds. 


5.2. QH - 

Our detection of C 4 H - constitutes the smallest reported col- 
umn density for this molecule in any source to date. Consistent 
with the trend described above for CgH - , the C 4 H - anion-to- 
neutral ratio of (1 .2 ± 0.4) x 10 -5 in TMC-1 (CP) is about nine 
times smaller than observed in L1527 by Agundez et al. (2008). 
Theory regarding C 4 H - chemistry is less well understood than 
for CgH - , but it seems plausible to again ascribe the lower value 
in TMC-1 (and the moderately low value of ( 8.8 ± 5.3) x 10 -5 
in Lupus- 1 A; Sakai et al. 2010), to chemical effects resulting 
from the lower densities and younger evolutionary states of these 
sources compared with L1527. 

The very low observed C 4 H - anion-to-neutral ratios com- 
pared with CgH - are at variance with chemical models that 
consider the formation of C 4 H - to be by radiative electron at- 
tachment to C 4 H at the rate calculated by Herbst & Osamura 
(2008). Recent ab initio calculations by V. Kokoouline et al. 
(2012, private communication) show that the C 4 H radiative 
attachment rate may be several orders of magnitude less 
than previously thought, and based on observations of the 
C 4 H - /C 4 H ratio in L1527, Agundez et al. (2008) calculated 
a C 4 H radiative attachment rate a factor of 122 less than the 
theoretical value of Herbst & Osamura (2008). However, using 
the laboratory data of Eichelberger et al. (2007) and Cordiner 
et al. ( 2012 ) identified that a significant pathway to smaller hy- 
drocarbon anions is via reactions of larger anions with atomic 
oxygen. Thus, an important source of C 4 H - is the reaction 

C 5 H - + O — > C 4 H - + CO. (1) 

As a consequence of this (and the analogous reaction that 
forms C 5 H - from CgH - + O), using the model of Cordiner 
et al. (2012) at «h 2 = 10 4 cm -3 (applicable to TMC-1), it is 
possible to set the radiative attachment rates for both C 4 H and 
C 5 H to zero and still produce an absolute C 4 H - abundance and 
C 4 H - /C 4 H ratio that are within a factor of a few of the observed 
values. It is therefore plausible that radiative electron attachment 
is not the dominant route to the formation of C 4 H - in molecular 
clouds. 


6 . CONCLUSION 

We conducted a search for the carbon-chain anion CgH - in 
seven nearby molecular clouds selected based on their strong 
HC 3 N emission. CgH - and CgH were detected in all sources, 
with a mean anion-to-neutral ratio of 3.0% (standard deviation 
0.92%). Combined with the four previous CgH - detections in 
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other sources, we deduce that anions are ubiquitous in the ISM 
wherever sufficient Q, H is present to be easily detectable. 

From the combined sample of 1 1 known interstellar sources 
with Q,FF detections, we confirm the trend identified by 
Cordiner et al. (2011) for anion-to-neutral ratios to be smaller 
in younger, less dense environments and larger in older, denser 
environments. This is consistent with the theory of Cordiner 
et al. ( 2012 ) that anion-to-neutral ratios are sensitive to the 
depletion of atomic O and H, and to the electron density, both 
of which become larger in denser, more chemically evolved 
cores. Molecular hydrogen number densities were derived for 
each source based on radiative transfer fits to observed HC 3 N 
spectra; a positive correlation is observed between Cf,H /Cf,H 
and «h 2 . 

We report the first detection of C 4 H in TMC- 1 , and derive 
an anion-to-neutral ratio C 4 H - /C 4 H = (1.2 ± 0.4) x 10 _5 (= 
0.0012% ± 0.0004%). This is the smallest anion-to-neutral ratio 
yet observed and confirms a problem with the theoretical value 
for the C 4 H electron attachment rate, which is apparently several 
orders of magnitude too large. The observed C 4 FR may instead 
be formed by reactions of the larger hydrocarbon anions ( C%, H 
and C 5 FC) with atomic oxygen. 

Given the widespread presence of anions in space, a com- 
bination of further dedicated laboratory, observational and the- 
oretical studies will be required in order to better understand 
their behavior and thus their full importance in astrochemistry. 
Measurements of radiative electron attachment rates for hydro- 
carbons, and branching ratios for reactions of anions with atomic 
oxygen will be of particular importance in this regard. 
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